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Abstract 

Spin degrees of freedom of charged nitrogen- vacancy (NV~) centers in diamond have large de- 
coherence times even at room temperature, can be initialized and read out using optical fields, 
and are therefore a promising candidate for solid state qubits. Recently, quantum manipulations 
of NV~- centers using RF fields were experimentally realized. In this paper we show; first, that 
such operations can be controlled by varying the frequency of the signal, instead of its amplitude, 
and NV~- centers can be selectively addressed even with spacially uniform RF signals; second, 
that when several NV~- centers are placed in an off-resonance optical cavity, a similar application 
of classical optical fields provides a controlled coupling and enables a universal two-qubit gate 
(CPHASE). RF and optical control together promise a scalable quantum computing architecture. 
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I. INTRODUCTION 



Impurity spins in diamond are among the most promising candidates for solid state 
quantum hardware. The so called (negatively charged) nitrogen- vacancy centers (NV~) 
have a low-lying spin triplet state 3 A with a large decoherence time (up to ~ 350 yus) at room 
temperature, which can be initialized and read out using a strong, spin- conserving optical 
transition to the excited 3 E state 01. The coherent manipulation of the 3 A state 

and its coupling to spins of 13 C [2], y] and N [4, 5] demonstrated the feasibility of NV~-based 
quantum devices. Though the direct coupling of different NV~- centers, necessary for a 
scalable architecture, would require placing them too close to each other (within a few nm), 

[inn 

coupling through an optical mode is possible [7|, |8|, [9j using Stark shifts, in order to tune the 
coupling on and off. (Stark shifts in NV~ were observed in bulk response 1(], 11, I2I as well 
as from individual centers [llj.) 

The use of local time-dependent fields for selective control is a natural approach, but it 
is not always easily achieved in the case of microscopic qubits. Here we suggest an approach 
which would allow us to address specific NV~- centers by tuning to their resonant frequency, 
which can be made position-dependent by the application of a static nonuniform magnetic 
field. We will also show that a similar approach using classical optical fields allows controlled 
coupling and universal two-qubit gates for NV~- centers in an optical cavity. 



II. MODEL 



An NV~- center is a negatively charged complex of a nitrogen impurity and a neighbouring 
vacancy. It can be formed as a result of nitrogen implantation in the diamond matrix; in 
experiments so far, the conversion from N to NV~ was achieved with a limited efficiency of 
about 5% Q]. It is therefore common to find an NV~- center close to a nitrogen impurity. 
Unlike an NV~- center, a nitrogen impurity does not have an electric dipolar moment and 



Let us choose 



does not couple to optical fields. Consider such a N-NV~ complex 
the [lll]-direction as the z-axis. The magnetic moment (spin 1) of the NV~- center and 
the eigenvalues of its z-component are S, M z = 0, ±1; for the N impurity (spin 1/2) they 
are denoted respectively by (l/2)a, m z = ±1/2; o is the vector of Pauli matrices. The 3 A 
ground state of the NV~-center is split by the crystal field, while the (M z = ±l)-states are 
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degenerate, and the (M z = 0)-state becoming the true ground state. The (M z = 0)-state 
is also the state leading to enhanced photoluminescence (through the excitation to 3 E and 
subsequent decay through a metastable level l A (see e.g. |l4|), allowing an optical readout 
{2, 3), 0, 0]. The external magnetic field along the z-axis splits the (M z = ±l)-states as well 
as the (m z = ±l/2)-spin states of the nitrogen impurity. The Hamiltonian of the system is 
(in the absence of an electric field) 

H = H NV + H N + H int , (1) 
where 

H nv = D(S z ) 2 + kB-S, (2) 
H N = -kB -a + Aa-I. (3) 



Here D = 2.88 GHz [l5j, ll6|, IjJ , k = 2.8 MHz/Gs the hyperfine splitting A = 86 MHz 



or 114 MHz depending on the position of the nitrogen in the lattice 18j, and / is its nuclear 
spin (1 = 1). The magnetic dipolar interaction 

H- mt = j[§-a-30-ft)(ff-fi)] (4) 

has a scale 7 ~ 6.5 MHz for the distance between NV~ and N of 2 nm; n is the unit vector 
in the direction connecting N and NV~. 

When B z = B ies = 514 Gs, the transition (M z = 0) ^ (M z = -1) in the NV _ - center 
is in resonance with the transition (m z = +1/2) <-> [m z = —1/2) in N. The term (j4j) then 
induces coherent transitions in the system, which were experimentally observed in 0, 0]. 
Other resonances, shifted by ~ 15 Gs to either side due to the hyperfine interaction in N, 
were also observed. (The hyperfine splitting in the NV~- center was too small to be resolved 
0]; we neglect it here.) 

In order to distinguish different NV~- centers, we now consider the application of a B z - 
field gradient. For example, if they are placed 10 /zm apart, a field gradient of 1 T/cm will 
produce a 30 MHz difference in the [M z = 0) <-» (M z = —1) transition frequency between 
the neighbouring centers, which is enough for our purposes, as we shall see below. 

III. RF CONTROL OF SINGLE-SPIN ROTATIONS IN NV - CENTERS 

If the field B z 7^ -B re s> the transitions (M z = 0;m z = 1/2) <-> [M z = —l;m z = —1/2) are 
suppressed. (Note that these energy conserving transitions do not conserve spin.) Single- 



qubit operation on a NV~- center can be then performed by applying a spatially uniform 
AC field along the y-axis, with the resonance frequency u y = D — kB z (of order of 1.5 GHz) 
corresponding to the (M z — 0) «-> (M 2 = — l)-transitions: 

Hy = liBy COSUytSy. (5) 

Due to the _B 2 -gradient, this frequency is different for different NV~- centers, and we have 
a frequency-based control. First, we go to the interaction representation: 

H y — > U NV H y U NV , (6) 

where 

Unv = exp{iH]yvt\ = (1 — S^) + (5^ cos nB z t + ifi^ sin nB z t) exp [zDt]. (7) 

After dropping the fast rotating terms (rotating wave approximation, RWA), we obtain the 
effective Hamiltonian 



tt _ kB v 



Sy [5*2, 5*; 



yl 



V2 



kB, 



y nv 



(8) 
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The operator in brackets acts as the Pauli matrix a y on the subspace {M z = 0, M z = — 1}. 
The AC field produces relatively fast rotations of a chosen NV~- center, in excess of 1 
MHz/G. 

The use of frequency instead of amplitude RF control is not dictated by the fact that 
the latter would require a local (within few microns) application of RF fields. The latter is 
feasible jljj] and, on a little larger scale, is being done in experiments with superconducting 
lux and phase qubits on a regular basis (see e.g. Q). Nevertheless the frequency control 



2l| has certain advantages including, but not limited to, less complex circuitry. 



IV. RF CONTROL OF NV -N COUPLING 

A relatively weak AC field with appropriate frequency can turn on the transitions 
(M z = 0; m z = 1/2) <-> (M z = —l;m z = —1/2). To see this, we again perform a uni- 
tary transformation of the Hamiltonian to the interaction representation, 



H - U N U NV HU ] NV U ] N - iU N U NV - 



d r 



UlyUU , (9) 
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where Unv was defined in ([7]), and 



U N = exp[iH N t] = cos [(kB z + AI z )t/2] + ia z sin [(kB z + AI g )t/2] 



(10) 



The resonance condition for the transition [M z = 0;m z = 1/2) «-> (M 2 = —l;m z = —1/2) 
is 2kB z + AZj — -D = 0. Assuming a detuning Su from resonance, i.e. 



2k5 2 + AI Z - D = 5co, 



we find (from Eq. (jUj)) the effective interaction 

rr /-i Q 2\ o I C0S(5^) 2 2 

# eff = 7(1 - 3nJS^ + 2 ^ 7(2 - 3n x - 3n y ) 



NV 



a,. 



h(t) 



(11) 



(12) 



(with the same notation as in Eq.®). Here h(t) denotes the rest of the terms, which all are 
fast rotating and should be dropped in the RWA. So is the second term in the r.h.s. of (1121) . 
unless the detuning Sou can be compensated. This is done by an additional field along the 
z-axis, B' z (t) = t]kB z sin out. (The corresponding term in the Hamiltonian is not affected by 
the transformation (J9])). After one more unitary transformation, 



at 



with 
U = exp 



-it] k B Z (S Z + ]-a z ) cos(ujt)/uj 



and assuming uo = 5u, we obtain, in the RWA, the following Hamiltonian: 

Ji (2t]kB z /lj 



H eS = 7(1 - 3n z )S z a z 



2V2 



-7(2 — 3n x — 3n 



^NV „ 1 „NV„ 

a x a x + a y a y 



(13) 



(14) 



(15) 



Otherwise the effective coupling is zero. Here J% is a Bessel function. This result, frequency- 
controlled NV~-N coupling, is confirmed numerically (Figs. [T|2|) . 

The coupling strength in resonance, 7 ~ 6.5 MHz, determines how far we should be 
off resonance in order to achieve decoupling. This frequency is large enough to justify the 
approximations leading to f|T5|) . The coupling is switched on by an AC field B' z {t). The 
coupling strength will be smaller than 7, but not drastically. Even choosing 5uo = 0.2 GHz, 
77 = B' z /B Tes = 0.05, and remembering that htB Tes w 1.5 GHz, we find the attenuation of the 
coupling strength w 0.32. It will still produce coherent transitions at ~ 2 MHz, which is 
fast compared to the decoherence times at room temperature as high as 0.35 ms. 
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FIG. 1: Maximum occupation probability for the NV - - center ground state, max , when 
initially in the excited state, as a function of the AC field frequency, u, for the NV~-N detuning 
5uj = 0.2 GHz and the AC field amplitude kB' z = 67 MHz (corresponding to rj = 0.05, with static 
field amplitude kB z = 1.34 GHz); 7 = 6.5 MHz. Note additional peaks at integral fractions of 6u, 
which correspond to the higher-order (multi-photon) processes, dropped from Eq. (|15p . 
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FIG. 2: (Color online) The occupation probabilities of the NV - - center ground state M z = 
(dashed line in blue) and the exited state M z = — 1 (in green) as a function of time, for the same 
choice of parameters as in Fig. [lj and u) = 0.2 GHz. The transition frequency (~ 2 MHz) agrees 
with the RWA value following from Eq. (|15p . 



As an aside, in the same way, one can show that when the system is in resonance, an 
AC field B' z (t) will suppress the transitions and freeze the spins in NV" and N (coherent 
destrustion of tunneling, see e.g. 22|]). This is not useful to the coupling scheme we are 
describing here. 

In order to perform arbitrary one-qubit rotations of nitrogen spins, as well as two-qubit 
gates between NV~ and N, we only need to initialize NV~ and perform single-qubit rotations 
on it (see e.g. [23(, where the phase qubit is analogous to an NV~- center, and the quantum 
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two- level system to a N impurity spin), which can be done using the RF control (Eq. (jSj)). 



V. NV -NV AND INDIRECT N-N COUPLINGS 



The scalability of the design requires the coupling between different NV - centers or 
NV~-N complexes. For macroscopic qubits this is done through the magnetic flux or charge 



coupling to the cavity modes 



23 



24 



25 



261 ] . In our case, unfortunately, the magnetic 



coupling is way too weak. Instead, we can use an optical cavity mode and two classical laser 



fields, along the lines of Refs. 



271 . This has the disadvantage of involving the 3 E state, 



where the decoherence rate is higher. On the other hand, the laser fields are easier to apply 
locally. By tuning the frequency of the laser field, the interaction strength can be controlled. 

Consider two NV~- centers placed in an off-resonance optical cavity. The Hamiltonian 
of the system is 



H = H Q + -ff add + H, 



L cavity j 



(16) 



H = lo c a) a +|r 



(/ + of) + (Zfil + ^7f ) + f (/ + erf ) + (*±^/ + 



2 '2 



(17) 



#fteld = #field,l + #field,2 = Qlx^X COs(u K " } t)(jf + g[ x ' 3 £[ COs(u\ L) t)7f 



(0), 



,(°) 



(-iM-i) 



x 
2 i 



(18) 



-"cavity -"cavity,! -^cavity, 2 



Here the operators 



g { iJe(a^ + a)af + g[ z L> e{a^ + a)r{ 



+fl£/e(a' +a)cr 1 + g\ z 'e(a' + a)r. 



(19) 



a] = \E,)(E,\ - \(A;M Z = 0) s )((A;M z = 0)0,|, 
if = I^X^-I - \{A-M Z = -l) j )((A;M z = -1),|, 
a* = \Ej)((A;M z = 0)A + \(A;M Z = 0),)(^-|, 
rj = |^-)((A;M 2 = + \(A;M Z = -1) 3 ){E 3 \ 

account for the optical transition 3 A <-> 3 E in the jth NV~- center. The term i/g e id de- 
scribes the interaction of the electric moments of the NV~- centers with the classical laser 
fields, S[°2 X \ with frequencies uf^ l \ polarized in the ^-direction; i/ cav ity describes their 



interaction with the optical cavity mode polarized in the z-direction; e is the "electric field 
amplitude for one photon in the cavity" . 

The idea of the approach remains the same as in the case of RF control of the NV _ -N 
coupling. For example, in order to induce the transition (A; M z = 0) <-> (E) in the NV~- 
center 1, we switch on the laser field Sf\ which is tuned to the frequency uj^ = oo c —E. After 
performing the unitary transformation with U = U caiV i ty Uo, where ?7 ca vity = exp[i / dt H caiVity ] , 
and Uq = exp[iH t], and then a RWA, the resulting term in the Hamiltonian will be, 

(0) (0)p(0) 

= -9 { Z(a^ + aaf), g% = , (20) 

n 

and similarly for the rest of the transitions. As in [27|, the coupling strength is proportional 
to the classical field amplitude. To target only one NV~- center, we now have two strategies. 
One is to reproduce our earlier approach and apply a non-uniform electric field. Then, due to 
the Stark shift, the resonance frequency of a given NV~- center will depend on its location, 
and the control is realized by applying uniform optical fields at specific frequencies. (Of 
course, due to the non-uniform static magnetic field applied to the system the (A «-> E) 
transition frequencies will already differ for different NV~- centers, but the difference is 
negligible.) 

The other strategy is to apply laser fields locally. Given the transition wavelength 637 
nm this may either put a lower limit on the spacing between NV~- centers, or require e.g. 
using evanescent modes in waveguides. 

The interactions « can produce single-^it rotations flQ. In onr S che m e they can 
be done more easily with RF pulses. On the other hand, two-qubit gates for different NV~- 
centers require long-range coupling, which can be achieved through virtual excitations in the 
cavity. For example, if we only apply the fields S[ and E% , and eliminate the cavity modes 
in linear order by the Schrieffer- Wolff transformation 28( , we obtain the effective interaction 

9 (0) (0) 

,^(00) / fl , eff,lfl'eff,2 / - + , - +\ / Q1 \ 

#eff,12 = —^( a l °2 + °2 H). (21) 

We we will also need transitions between the 3 A levels M z = 0(— 1) and the 3 E state. They 
can be realized by applying another laser field with rc-polarization and at a corresponding 
resonant frequency ~ 470 THz (A = 637 nm). The cavity degrees of freedom or other NV~- 
centers will not be involved. 
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Using (12T1) . we can produce, e.g., an entangling transformation 
(a\M z = Q)!+/3\M Z = -l)i)®|M, = 0) 2 -> (a|M 2 = 0) l( g)|M 2 = -1) 2 +/?|M 2 = -1>!®|M, = 0) 2 ).(22) 
To achieve this, we perform the following set of operations: 

1) 7r-pulse between \M Z = 0) 1 and \E) 1 ; the result is (a\E} l + (3\M Z = ® |M 2 = 0) 2 ; 

2) vr-pulse of the interaction ([21]), resulting in (a\M z = 0) 1 ®\E) 2 +(3\M Z = -l) 1 ®\M z = 0) 2 ); 

3) vr-pulse between \E) 2 and \M Z = -1) 2 ; the outcome is (a|M a = 0)i<g>|M 2 = -1} 2 +/3\M Z = 
— l)i (g> |M 2 = 0)2), as required. After enabling the effective NV~-NV~ coupling through the 
cavity, the operations on N impurities coupled to different NV~- centers can be realized in 



the same way as in Ref. 23 1. 



Although the above procedure nicely demonstrates the possibility of using the cavity to 
perform two-qubit gates, it takes any state of the form (a\M z = 0)i+@\M z = — l)i) <S> \M Z = 
—1)2 outside the computational basis. Instead, the C-phase gate could be implemented as 
follows: 

1) 7r-pulse between \M Z = 0) 1 and \E)i, 

2) 27r-pulse of the interaction fT2"Tl) ; 

3) 7r-pulse between \M Z = 0)i and \E)\. 

An inspection of the above procedure shows that the three states \M Z = 0)i|M 2 = 0) 2 , 
\M Z = —l)i\M z = 0)2 and \M Z = —l)i\M z = — 1) 2 are left unchanged at the end of the 
procedure, whereas the state \M Z = 0)i\M z = — 1) 2 acquires a minus sign (note that we are 
not including here the phases accumulated as a result of single-qubit Larmor precession). 
This two-qubit gate, along with single-qubit rotations, form a universal set of gates for 
quantum computing. 

The requirements to the optical cavity are high, but not impossible. In order to resolve 
~ 1.5 GHz against the ~ 470 THz resonant frequency, the quality factor of the cavity should 
be at least of order 3 x 10 5 . We could somewhat improve the situation by adding one more 
step and swapping the states M z = and M z = 1. This can be done again using the RF 
field B y (t) (Eq. (jSJ)), this time with the frequency u' = D + kB z (of order of 4.5 GHz). This 
increases the difference of the states involved in the optical coupling from ~ 1.5 GHz to ~ 3 
GHz. 
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VI. CONCLUSIONS 



We propose a frequency-controlled approach to coherent manipulation of spin states of 
NV~- centers and NV~-N complexes in diamond. It allows to address different spins through 
the difference in their resonance frequencies, induced by a static non-uniform magnetic field. 
The time-domain manipulations are performed using uniform RF fields. Different NV~- 
centers and NV~-N complexes can be coupled optically through the virtual excitations in an 
optical cavity. Here both frequency control with spatially uniform AC fields and with local 
AC fields are possible. The required cavity quality factor is high, but achievable. Our results 
show that small-scale quantum information processing devices based on impurity spins in 
diamond may be feasible in the near future. 
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